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\  I  Abstract 

ihis  paper  deals  with  a  nonparametric  subset  selection  procedure 
for  a  two-way  layout  problem.  The  treatment  effect  with  the  largest 
unknown  value  is  of  interest  to  us.  The  block  effect  is  a  nuisance 
parameter  in  this  problem.  The  proposed  procedure  is  based  on  the 
Hodges- Lehmann  estimators  of  location  parameters.  The  asymptotic  relative 
efficiency  of  the  proposed  procedure  with  the  normal  means  procedure  is 
evaluated.  It  is  shown  that  the  proposed  procedure  has  a  high  efficiency. 


Key  words:  Nonparametric  procedure,  subset  selection,  two-way  layout, 
treatment  effect,  nuisance  parameter,  asymptotic  relative  efficiency. 
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NONPARAMETRIC  SELECTION  PROCEDURES  FOR 
A  TWO-WAY  LAYOUT  PROBLEM* 

by 

Shanti  S.  Gupta  and  Lil-Yuh  Leu 
Purdue  University 


1.  Introduction 

Consider  a  two-factor  complete  block  design  with  one  observation  per 
cell.  Let  the  observable  random  variables  be  X^,  i  =  l,2,...,k;  a  =  1,2,... 
and  consider  the  linear  model 


(1.1)  e1a.  ^  ef  •  0. 

where  X1a  Is  the  observation  under  treatment  i  in  the  ath  block,  u  is  the 
mean-effect,  is  the  effect  of  treatment  i,  Ba  is  the  block  effect  for 

the  ath  block  (nuisance  parameter),  and  the  e^a,  a  =  1,2 . n  are  error 

components.  It  Is  assumed  that  the  error  components  are  Independent  and 
Identically  distributed  with  a  continuous  cumulative  distribution  function 

L 

(cdf)  F(e) ,  e  €  1R  (the  real  k-space),  where  F(e)  is  symmetric  in  its 

L 

arguments,  that  is, for  any  e  €  F  and  any  permutation  (11,...,ik)  of 
(1,...,k),  we  have 


(1.2) 


F(e, «... ,et)  ■  F(e 


Ej  ). 


4  »  •  •  • 

'1  *k 


Let  <...<  be  the  ordered  e^'s.  Suppose  that  we  are  interested 
in  the  treatment  with  the  largest  unknown  parameter  (if  more  treatments 
than  one  have  equal  to  then  exactly  one  of  these  treatments  is 
"tagged1*  as  the  best  treatment).  Correct  selection  denotes  the  selection  of 
aqy  subset  containing  the  population  with  (or  the  "tagged"  population). 


*this  research  was  supported  by  the  Office  of  Naval  Research  Contract 
N00014-75-C-0455  at  Purdue  University.  Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of  the  United  States  Government. 


For  the  nonparametric  approach  to  the  one-way  layout  problems,  most  previous 
authors  have  considered  procedures  based  on  a  class  of  rank  order  statistics. 

It  was  pointed  out  by  Rizvi  and  Woodworth  (1970)  that  there  are  difficulties 
associated  with  these  procedures,  mainly  because  the  least  favorable  configura¬ 
tion  is  usually  not  known.  Randles  (1970)  and  Ghosh  (1973)  have  considered 
the  procedures  for  the  one-way  layout  based  on  the  Hodges-Lehmann  estimators. 

It  was  shown  by  them  that  the  procedures  based  on  the  Hodges-Lehmann  estimators 
have  high  efficiency.  In  this  paper  we  use  the  results  of  Puri  and  Sen 
(1967)  to  derive  a  subset  selection  procedure  for  the  largest  unknown 
parameter 

2.  Robust  Compatible  Estimation 

In  the  model  (1.1),  for  1  <  i  <  j  <  k,  let  X..  =  X.  -  X.  , 

—  *“*■  IJ»0L  lot  Jut 

e.4  =  e.  -  e.  ,  a  =  1,2 . n  and  =  e4-e.,  then  for  a  fixed  a, 

ij,a  la  Ja  1j  1  J 

the  ath  block,  we  can  write 

(2.1)  X.  i  *  a.  i  ^  e . .  . 

'  1j»a  ij  lj,a 


From  Assumption  (1.2),  e,4  have  common  distribution,  say  G,  which  is 

l  j  ,a 

symmetric  about  zero.  Hence  X.j  j,...,X.j  n  are  i.i.d.  with  common 

cdf  G(x-Aij),  for  1  ^  1  <  j  <_  k.  We  assume  that  G  is  continuous,  but 

otherwise  unknown.  Let  R.,4  *  Rank  of  |X<4  I  among 

i  j  ,a  ij»a' 

lxij,1|,..,lx1j,n!  and  let  X^  *  (X^j . xij,n^  Cons1der  the  one- 

sample  signed  rank  statistic 

<2-2>  h14.n«?1j>  ’  n‘1  bnjn,« 

a  i 

Where  Z„  Is  either  one  or  zero  as  follows:  If  the  ath  smallest  observation 
n,a 

among  |X^j  j  |,...,|X^n|  .corresponds  to  X^t  (for  some  t),  then 

1  «  1  if  X44  t  >  OorO  if  X14  t  <  0.  E_  „  Is  the  expected  value  of  the 

n,a  ij.t  1J,t  n,a 


ath  order  statistic  of  a  sample  of  size  n  from  a  distribution  4>*(x)  given 
by 


(2.3) 


**(x) 


I 


rl>(x)  -  *(-x) 
0 


If  x  >  0 


If  x  <  0. 


Throughout  this  paper,  we  shall  assume  that  tj*(x)  and  G(x)  satisfy  the 
following  assumptions  (see  Purl  (1964)  and  Purl  and  Sen  (1967)): 

(1)  *(x)  is  symmetric  about  x  ■  0,  that  Is  <j»(x)  +  1<(-x)  ■  1. 

<2>  w  j,[Er.«-*',(ST>’Zn.a- Vn_i)- 

(3)  J(u)  *  ^(u),  0  <  u  <  1  Is  absolutely  continuous  and 

|jW(u)|  -  |d1J(u)/du1 1  <  HCu(l-u)]”1^  +\  1  -  0,1,2,  for  some 
M  and  some  a  >  0. 

(4)  G  Is  a  continuous  cdf,  differentiable  In  each  of  the  open  intervals 
(-•.tj);  (a^ .ag) . .  «(*s-i »as)  *  (®5*"°^*  some  a-|,...,a^  and  the 
derivative  of  G  Is  bounded  In  each  of  these  Intervals. 


(5)  The  function  ^  J(G(x) )  Is  bounded  as  x  ±». 


It  Is  easy  to  see  that  hi j ,n(x1J  l+a,,,*,x1j,n+a^  1s  a  non"decreasin9 
function  of  a  for  fixed  x^j  and  when  ■  0,  the  distribution  of  h-jj>n  is 
symmetric  about  a  fixed  point  u  *  y  EjV|,  where  V  has  cdf  <j>. 

Let 


**1i 

■  supU: 

C2.4) 

#* 

■  InfU: 

and  let 

*  1 
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Now  from  Hodges  and  Lehmann  (1963),  is  a  translation  invariant  robust 
estimator  of  a^j  and  has  a  distribution  which  is  symmetric  about  a^.  Note 
that  the  estimates  a^  are  Incompatible  (see  Lehmann  (1964))  in  the  sense 
that  they  do  not  satisfy  the  linear  relations  satisfied  by  the  differences 
they  estimate.  This  leads  to  certain  ambiguities.  To  derive  compatible 


estimators,  let 


(2.5) 


-  1  k 

Ai  =  TT  J,  *1i#  *1i  H  0  for  1  *  1,2 . k* 


’j 

Then  by  minimizing  T  (L.-a..)  with  respect  to  e's,  we  obtain  the 
ifj  1J 

compatible  or  adjusted  estimators  of  a^  as 

(2.6)  Z^.  *  A^~Aj ,  i  t  j. 

k 

Note  that  E(a^)  ■  0^  since  J  *  0,  hence  E(Z^ ^ )  *  a.j.  Puri  and  Sen 
(1967)  have  proved  the  following  theorem: 

Theorem  2.1.  The  joint  distribution  of  {n^fZ^-A  ik);  1  *  1  ,2 . k-1)  is 

asymptotically  normal  with  zero  means  and  a  covariance  matrix 
r  *  (Yij).i.J  *  1.2 . k-1  where 


(2.7) 


2oq  If  1  -  j 


4  if  MJ 


and  oq  »  [A2  +  (k-2)Xj(G)]/kB2, 


where 


(2.8) 


j  J2(u)du, 


B  -  /  fy  J(6(x))dG(x), 


and  Xj(G)  =  /  /  J(G(x))J(G(y))dG*(x,y), 


5 


G*(x,y)  Is  the  joint  cdf  of  e^j  a  and  ei£  a(j  f  t)  whose  marginal  cdf's 
are  G(x)  and  G(y),  resect Ively. 

Moreover,  using  the  translation  Invariant  property  of  we  have 

I  J 

the  following  lemma: 

Lemma  2.2. 

(1)  A1(x11-c1,...,xln-c1,...,xkl-c|( . ^n-Cfc) 

=  Vxll . xln***‘  *xkV*  ,xkn^  +  c_ci* 

-  1  k 

where  c  *  £  I  c4. 

*  1*1  1 

(2)  The  distribution  of  Is  independent  of  e. 

3.  A  Nonparametric  Procedure  for  Selecting  the  Best  Treatment 

Based  on  the  estimators  defined  in  (2.5),  we  propose  a  selection 
procedure  as  follows: 

(3.1)  R,:  Select  treatment  1  iff  L  >  max  A.-d,, 

'  1  ~  l<j<k  J  1 

where  dj  >  0  Is  determined  so  as  to  satisfy  the  basic  probability  requirement. 

Let  n  *  (e  *  (e1,...,ek)}  be  the  parameter  space  and  CS  stand  for  a 

correct  selection  which  means  that  the  selected  subset  contains  the  best 

treatment.  For  a  given  constant  P^k’1  <  P*  <  1),  the  basic  probability 

requirement  Is 

Inf  Pfl(CS|R,)  >  P*. 

e€n  2 

Let  Aj-^  <...<  Aj*kj  denote  the  ordered  a^'s  and  a^  denote  the 
unknown  estimator  associated  with  the  parameter  1  <  1  <  k.  If 
Pj(8|Rj)  denotes  the  probability  that  the  treatment  (j)  is  selected 
(treatment  (j)  Is  associated  with  parameter  e^).  We  have  the  following 


Lemma  3.1.  Pj(0|R-|)  is  increasing  in  when  all  other  components  of 
9  are  fixed,  and  is  decreasing  in  i  j  j,  when  all  other  components 
of  e  are  fixed. 


Proof.  PjtjlR,)  ■  P,(i0)  >  ““(D'O,) 

=  P|(i(i)‘6[1]'i(j)+e[j]  -  dl  *  e[j]‘e[1]>  '  +  1  "  ' . k)- 

By  Lemma  2.2,  the  distribution  of  i )“e[i 3~^( j )  +  6[j]»  1  4  3»  1  *  l,2,...,k 
is  independent  of  e.  Hence  Pj(e|R^)  is  increasing  in  0^  an<l  Is  decreasing 
in  0[i]»  i  t  J*  Note  the  above  property  is  usually  called  strong  monotonicity 
property. 

Corollary  3.2.  inf  P  (CS|R,)  *  inf  P  (CS|R,)  which  is  independent  of 

/» ^  O  I  n  « n  W  | 


0€R  - 


=...=  9k  =  0,  where  nQ  =  {e  €  nje^  *...=  e^}, 


Proof.  The  proof  follows  from  Lemma  3.1. 

For  large  sample  we  can  define  d1  as  in  Theorem  3.3  given  below. 

Theorem  3.3.  For  given  P*(k  <  P*  <  1).  If  Oq  <  we  have 

(3.2)  d-j(n)  =  n"^  daQ  +  o(n"^)  as  n  +  », 
where  d  Is  the  solution  to  the  equation 

(3.3)  Q(d/ /?,..., d/v^  »  P*, 

Q  is  the  joint  cdf  of  a  normally  distributed  vector  (V1,...,Vk_1)  with 

(3.4)  EfVj)  -  0,  VarfV^  *  1  and  Cov(¥1,Vj)  =  1/2,  1  *  j. 


7 


Proof.  11m  Inf  P0(CS |R^)  *  11m  pq(^(i )”^(k)  -  di^n^»  1  =  1,2,. ...k-l) 
ih®  ^  rn® 

a  P(V*  £  11m  r£d.(n)//?  Oq,  1  =  1,2 . k-l),  by  Theorem  2.1, 

n-»® 

where  Pq  denotes  the  probability  Is  computed  under  =...=  ek.  If  d 
Is  the  solution  of  (3.3),  then 

d-j(n)  =  n"^do0  +  °(n"^)  as  n  -*  ». 

00  k  l 

Remark:  The  solution  of  (3.3)  Is  also  a  solution  of  /  $  (x+d)  d$(x)  =  P*, 

•oo 

where  ♦  Is  the  cdf  of  standard  normal.  This  has  been  shown  by  many  authors 
(see  for  example,  Gupta  (1963)). 


Determination  of  the  Minimum  Common  Sample  Size 

Let  E(S | R-. )  denote  the  expected  size  of  the  selected  subset  using 
1  k 

rule  R-j,  then  E{S|Rj)  *  X^Pj(ejR-j).  Having  determined  d-j(n)  from  (3.2), 
one  may  determine  the  common  sample  size  n  by  Imposing  the  additional 
requirement  that  E(S|Rj)  £  1+e,  for  some  e  >  0,  whenever  e  lies  in  a  given 
proper  subset  of  a ,  for  example,  the  subset  defined  by 


(3.5)  epj  *...*  =  -  5*»  <5*  >  0* 

It  will  be  convenient  In  the  sequel  to  replace  (3.5),  when  the  sample  size 
Is  n,  by 

(3.6)  8[i]  ■•••"  e[k-1]  *  0[k]  "  s(n)- 
(See  Sartlett  and  Govlndarajulu  (1968)). 

Theorem  3.4.  For  given  e  >  0,  with  d-j(n)  given  by  (3.2)  and  n  determined 
by  E(S|R<|)  £  1+e  for  e  satisfying  (3.6).  Then  as  n  -*•  », 

(3.7)  -  •***€■'  /«0  +  o(n^). 
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where  c(e)  is  the  solution  to  the  equation 

(3.8)  Q((c+d)/»'?,...,(c+d)/v'2)  +  (k-l)Q(d/^ . d/^,(d-c)/^f)  =  1+e, 


where  Q  is  defined  as  in  Theorem  3.3. 


Proof.  E(S  1^)  =  l  Pe(r#(A(i)-A(j))/^a0  <  r^d^n)/^,  i  =  l,...,k-l, 

1  t  j.  ^(i(kri(j)-5<n))/^»0i  r^(<i,(n)-«(">)/<?»0) 

♦  Pe(r^(£(1)-4(|()+«*n*)/«'i5j0  <  i£(d1(n)«*n*)//S>0, 


If  e  satisfies  (3.6),  then 

lim  E(S|R1)  =  k'l  P(Vi  <  lim  n£d ^n)//^,  i  =  l,...,k-l,  i  +  j , 
rv*»  j=l  n-*® 

V.  ,  £  nF(d1(n)-6^n^)/i^a0)+P(Vi  <_  lim  ^(d^nJ+d^^/^bQ, 

i  =  1,2 . k-1) 

=  (k-l)Q(d/v'?,...,d/v'?,  (d-c)/v^)+Q((d+c)/^,...,(d+c)/v^)  *  1+e. 

Hence  c(e)  is  the  solution  of  (3.8)  iff 

6(n)  =  n~^c(e)og  +  o(n“^)  as  n  +  ®. 


Remark:  The  common  sample  size  n  required  to  satisfy 

Inf  P  (CS|R,)  =  P*  and  E(SjR,)  <_  1+e  for  e  satisfying  (3.5)  is 
6€0  - 

(c(e)o0/5*)2.  Note  that  n  is  a  function  of  k,  P*,  6*,  and  e. 


4.  A  Selection  Procedure  for  the  Normal  Case 

If  we  assume  that  (e1  ,...,ek  )  are  jointly  normally  distributed  with 

zero  means  and  the  covariance  matrix  o  (  *.  ),  a  *  1,2,. ...n  where 

p  *1 

2 

-l/(k-l)  <  p  <  1  is  known  and  o  <  »,  may  be  known  or  unknown.  The  usual 


9 


least 


1  ” 

square  estimator  of  e,.  would  be  X.-X  where  X.  =  ^  7  X. 

i  i  i  n  ia 

a=  I 


1  K 

X  *  ^  l  X...  It  is  easy  to  show  that  the  vector  (X-j-X,. . .  ,Xk~X)  has  a 
joint  normal  distribution  with  mean  vector  (6p...,ek)  and  the  covariance 


matrix 


2n  .  /1c  "IT 
<4-')  2jrai  i  i 


1  1 


and  hence  the  vector  (X-j  -  has  a  j°int  normal  distribution 

with  mean  vector  (e1-ek,. .. ,ekl-9k)  and  the  covariance  matrix 


(4.2) 


1 

1  **., 


& 


s 


We  assume  that  o c  is  known  and  propose  a  selection  procedure  R2  by 

(4.3)  R0:  Select  treatment  i  iff  X..  >  max  X.  -  d0. 

1  1  ”  l<j<k  J  1 

It  is  easy  to  see  that  inf  Pfl(CS|R?)  *  inf  P  (CS|R?)  and  is  independent 

9€fi  -  0€fio  ■  - 

of  *...*  ek  *  e.  Thus,  similar  to  Theorem  3.3  and  Theorem  3.4,  we 
have  the  following  theorems: 

Theorem  4.1.  For  given  P*  (k-1  <  P*  <  1),  if  we  let 

Inf  P  (CS|R«)  *  P*,  then 
e€fi  - 

(4.4)  dg(n)  =  n^di/Pp'  o, 
where  d  is  the  solution  of  (3.3). 

Theorem  4.2.  For  given  e  >  0  and  d2(n)  *  n^d^T^T  o,  If  n  is  determined  by 
E(S|Rg)  <,  l+t  for  e  satisfying  (3.6),  then  as  n  -►  ~ 


(4.5) 


<5^  =  n  8  c(e) /1-p  o, 

where  c(e)  Is  the  solution  of  (3.8). 

Suppose  that  the  joint  distribution  F  of  (ela,. .. ,eka)  is  unknown, 
but  the  variance  of  e-ja  is  finite.  By  the  central  limit  theorem,  the  joint 

distribution  of  {r#(X.-X  k~ei+ek^’  *  =  ^*2 . k-l>  is  asymptotically  normal 

with  zero  means  and  the  covariance  matrix  (4.2).  We  can  still  use  the 
procedure  R2  given  by  (4.3).  For  large  samples  we  have 
d2(n)  =  n  8  dA-p  a  +  o(n  8)  and  6^  =  n"^c(E)/1-p  a  +  o(n”^),  where  d  is 
the  solution  of  (3.3)  and  c(e)  is  the  solution  of  (3.8). 

For  any  two  procedures  R-|  and  R2  satisfying  the  basic  P*-condition,  let 
us  define  the  asymptotic  relative  efficiency,  say 

ARE(R, ,R9)  =  lim  nD  (e)/nD  (e)  for  the  given  parametric  configuration  (3.5), 

1  6  e+0  K2  K1 

where  nD  (e),  i  =  1,2  are  the  sample  sizes  required  to  achieve  the  same 
Ki 

expected  size,  1+e.  Then  we  have  the  following  theorem: 

Theorem  4.3.  ARE( R] ,R2)  *  {2o2(l-P)B2/A2{[kA2/2[A2+(k-2)xJ(G)]},where 
2 

A  ,  B  and  Aj(G)  are  defined  in  (2.8). 

Proof.  For  procedure  Rj»  putting  6^  =  6*,  from  (3.7)  we  have 
nR  (e)  •  (c(e)o0/5*)2, 

where  c(e)  is  the  solution  of  (3.8).  (Note  that  e  +  0,  then  n  <»). 
Similarly,  for  procedure  R2,  we  have 

nR  (e)  =  (c(e)/1-p  o/6*)2. 
c 

Hence 

ARE( R, ,R2)  *  (l-p)o2/a2 


11 


*  {2b2(l-p)B2/A2HkA2/2[A2+(k-2)Xj(G)]}. 

Remarks : 

(1)  Barlow  and  Gupta  (1969)  define  ARE(R-|,R2)  =  lim  nR  (e)/nR  (e)  for  the 

c  e+0  c  *1 

given  parametric  configuration  (3.5),  where  nD  (e),  i  =  1,2  are  the 

Ki 

sample  sizes  required  to  achieve  the  same  expected  size  (say  e)  of  non¬ 
best  populations  selected.  If  we  consider  the  case  where  expected  size 

refers  only  to  the  number  of  non-best  populations  in  the  selected  subset, 

2 

we  have  nR  (e)  =  (c'(e)o0/6*)  ,  where  c'(e)  is  the  solution  to  the 
equation 

(k-l)Q(d//?,. . . ,d/V*?,  (d-c')//?)  *  e. 

Similarly,  we  have  nD  (e)  =  (c,(e)/Pp’  o/6*)  ,  and  hence 

k2 

ARE(RrR2)  *  (l-p)o2/c2 
which  Is  the  same  as  In  Theorem  4.3. 

1  2 

(2)  Puri  and  Sen  (1967)  proved  that  Xj(G)  A  ,  hence 

kA2/2[A2+(k-2)Xj(G)]  >  1  and  ARE(R^ ,Rg)  >  2(l-p)o2B2/A2.  The  variance 
9  2  2  2 

of  G  Is  2(1-p)o  ,  hence  2(l-p)o  B  /A  is  the  ARE  of  the  one-sample  rank 
order  tests  (for  location)  with  respect  to  the  Student's  t-test  when  the 
parent  distribution  is  G(x).  If  we  use  the  normal  scores  estimator,  we 
have  ARE(Rj ,R?)  ^  1.  If  we  use  the  Wilcoxon  scores  estimator,  then  for 
any  F,  we  have  ARE(R^,R2)  >.  0.864  and  ARE(R-|,R2)  s  3/u  when  F  is  normal. 
Hence  the  procedure  given  by  (3.1)  has  "high"  efficiency. 

5.  Relative  Performance  of  R-|  and  Rg 

In  Section  4,  we  consider  the  parameter  points  satisfying  (3.6).  When 
the  condition  (3.6)  is  not  satisfied,  but  the  ratio  of  sample  sizes,  m  for  R^ 
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n  p  o 

and  n  for  R«,  satisfies  lim  =  (1-p)a  /of,  then,  for  large  n,  the  procedures 
c  n-*»  m  u 

and  R2  have  approximately  the  same  probability  of  a  correct  selection  and 
expected  size. 

n  p  p 

Theorem  5.1.  Let  n  and  m  =  g(n)  satisfy  lim  ^  =  (l-p)o  /og,  then  the  procedures 

n-*» 

and  R.j  have  the  same  asymptotic  probability  of  a  correct  selection  and  the 
same  expected  size  for  any  parametric  configuration. 

Proof.  For  procedure  R2,  consider  any  sequence  of  parameter  points  satisfying 

6[k]  ‘  e[i]  =  6in  =  n’2y^  a5i  +  °(n"2)» 
i  ■  l,...,k-l  and  for  some  i,j,  fi.  4  6*.  <s .  4  0-  Then 

'  J 

lim  P  /n\(CS|R2)  =  P(V.  <  lim  i^(dz(n)+5i  )//2'(1'-P)  o,  i  =  1 . k-1) 

n-x»>  0'  n->«> 

=  P(Vi  <  (d+«.)/v^,  i  =  1,2,..., k-1) 

-  Q((d+61)/^,...,(d+6k-1)/^) 
and  for  1  <_  j  <  k-1 , 

lim  P.(e(n)|R2)  =  P(Vi  1  (d+«r«j)/^,  i  =  1,2 . k,  i  f  j) 

n-x» 

■  Q((d+61-6j)/^,...,(d+6k-6j)/^). 

For  procedure  R^,  m^a0  -v  r\^A^  a,  so 

e[kl  ‘  6[i]  =  "'Vi  +  o(m^}  *  5im*  1  =  1,2 . k”1’ 

Hence 

Urn  P  /  %(CS|R-j)  *  P(Vi  <  (d+6^//?,  i  =  1,2 . k-1) 

nH«  e'  ' 

*  Q((d+«j )/»#,. . .  ,(d+6k _ -j ) / - 

and  for  1  <  j  <_  k-1. 
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llni  Pj(e^|Ri)  *  Q((d+«1-fi  .)//? . (d+6k-6  .)//?). 

RVtoo  j  J  J 


In  the  above  parameter  points,  we  assume  that  tend  to  zerc 

at  the  n”^  rate.  If  any  difference  tends  to  zero  more  rapidly,  we  replace 


by  0,  and  if  it  tends  to  zero  more  slowly,  or  tends  to  a  finite  limit, 
then  we  replace  6^  by  ®,  and  still  obtain  the  same  asymptotic  behavior.  This 
completes  the  proof  of  the  above  theorem. 


6.  Estimation  of  B  and  Xj(G) 

In  practical  application,  for  large  n,  the  procedure  Rj  can  be 
rewritten  as 

don 

(6.1)  R, :  Select  treatment  i  iff  a .  >  max  a.  - -, 

1  1  ~  l<j<k  J  ^ 

*  L  ] 

where  d  is  the  solution  of  (3.3)  or  /  *  (x+d)d«(x)  =  P*.  However,  aQ 

2 

is  still  unknown.  We  need  to  find  a  consistent  estimator  of  Oq.  Since 
oq  ■  {A2+(k-2)Xj(G)}/kB2,  where  A2  =  ^J2(u)du  is  known,  but 

B  ■  /  ^  J(G(x))dG(x)  and  Xj(G)  =  /  /  J(G(x))J(G(y))dG*(x,y)  are 

—09  —CO  —00 

unknown.  Our  problem  is  to  find  consistent  estimator  of  B  and  Xj(G), 

The  consistent  estimators  of  B  and  Xj(G)  can  be  found  from  (4.7) 

and  Theorem  4.2  of  Puri  and  Sen  (1967);  let  these  be  Bn  and  Ln,  respectively. 

A  o  o  a  *2  2 

Then  oQ  *  [A  +(k-2)Ln]/kBn  is  a  consistent  estimator  of  Oq.  Hence,  for 

large  n,  the  procedure  R^  is  defined  by 

«  doQ 

(6.2)  R.:  Select  treatment  1  iff  Aj  >  max  a4  -  — -. 

l<,j<k  ^  /n 

Remark;  If  *(x)  is  the  cdf  of  U(-l,l),  then  J(u)  *  2u-l,  0  <  u  <  1,  hence 
A2  ■  1/3.  If  G'(x)  ■  g(x)  exists,  then  B  *  2/g2(x)dx  and 
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Aj(G)  *  4/  /  G(x)G(y)dG*(x,y)-l.  It  has  been  shown  by  Doksum  (1967)  that 
-00  -00 

OO  CO 

/  /  G(x)G(y)dG*(x,y)  can  be  estimated  by  n'/n-TWn-glkfk-l  (number  of 


sixtuples  (1,j,*,a,e,y)  with  i,j,£  distinct;  a,B,y  distinct,  and 

VXja  <  Xi8-Xj8'  VX*a  <  ^y  -  V* 
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